Competing Orders and Spin-Density- Wave Instability in La(Oi_ x Fa;)FeAs 



J. Dong, H. J. Zhang, G. Xu, Z. Li, G. Li, W. Z. Hu, D. Wu, G. F. Chen, X. Dai, J. L. Luo, Z. FangQand N. L. Wan^l 

Beijing National Laboratory for Condensed Matter Physics, Institute of Physics, 
Chinese Academy of Sciences, Beijing 100080, People's Republic of China 



00 

o 

O 

(N 



(N 



C ■ 
O 1 

o , 

u ■ 

Oh 

s : 

a ■ 



-a 
c 

o 
o 



> 

(N 

cn 
rn 

o 

oo 
O 



X 



The interplay between different ordered phases, 
such as superconducting, charge or spin ordered 
phases, is of central interest in condensed matter 
physics. The very recent discovery of supercon- 
ductivity with a remarkable T c = 26 K in Fe-based 
oxypnictide La(Oi_ 2 ,F 2; )FeAs is a surprise to the 
scientific community [1]. The pure LaOFeAs it- 
self is not superconducting but shows an anomaly 
near 150 K in both resistivity and dc magnetic 
susceptibility. Here we provide combined exper- 
imental and theoretical evidences showing that 
the anomaly is caused by the spin-density-wave 
(SDW) instability, and electron-doping by F sup- 
presses the SDW instability and recovers the su- 
perconductivity. Therefore, the La(Oi_ 2 ,F a; )FeAs 
offers an exciting new system showing competing 
orders in layered compounds. 

Both charge-density wave (CDW) and spin-density- 
wave (SDW) instabilities may develop in the presence 
of Fermi surface (FS) nesting[2j, as shown in many tran- 
sition metal compounds, such as 2H-NbSe2[l] and CrQ. 
The difference is that CDW couples to lattice and SDW 
couples to spin. The CDW or SDW instabilities may 
compete with other possible orderings, and complicated 
phase diagrams are often drawn due to their interplay. 
This is exactly what happens for LaOFeAs as shown in 
this paper. LaOFeAs crystalizes in layered square lat- 
tice with Fe layers sandwiched by two As layers (up and 
down), each Fe is coordinated by As tetrahedron. Its 
electronic properties are dominated by the (FeAs)-triplc- 
layers, which contribute mostly to the electronic state 
around Fermi level (Ef). The electronic structure of 
LaOFeAs is quasi-two-dimensional and very similar to 
typical semi-metal, namely there are hole-like FS cylin- 
ders around the T-Z line of the Brillouin zone (BZ), 
and electron-like FS cylinders around the M-A line of 
the BZ (except a small 3D FS around the Z point of 
BZ) [H, [o, 0, @, 0] • We will show here that surprisingly 
strong FS nesting exists by connecting the hole and elec- 
tron FS by a commensurate vector q=(ir, w, 0). This 
leads to SDW instability, and it is the main cause of the 
anomaly observed at 150 K experimentally. The spon- 
taneously symmetry-broken SDW state is characterized 
in terms of reduced carrier density due to (partial) FS 
nesting, enhanced conductivity due to the reduction of 
scattering channel, and loss of 4-fold rotational symme- 
try with negligible change of lattice. 

Series of La(0i_ 2 ,F a; )FeAs samples are fabricated, and 
Fig. 1 (a) shows the temperature dependence of the re- 



sistivity. The pure LaOFeAs sample has rather high dc 
resistivity value and very weak temperature dependence 
at high temperature, but below roughly 150 K, the resis- 
tivity drops steeply, with an upturn at lower temperature 
(below 50 K). At 2% F-doping, the overall resistivity de- 
creases and the 150 K anomaly shifts to the lower tem- 
perature and becomes less pronounced. At 3% F-doping, 
the anomaly could not be seen, and a superconducting 
transition occurs at T c =17 K. With further increasing 
F content, and the superconducting transition tempera- 
ture increases with the highest T c =28 K seen at 6~8% 
F-doping. A T c vs. F-content phase diagram is plotted 
in the inset of the figure 1 (a) . 

Apparently, superconductivity competes with the 
phase showing the anomaly in this system. To identify 
whether the anomaly at 150 K in LaOFeAs is due to 
a phase transition, we performed the specific heat mea- 
surement. Figure 1 (b) shows the specific heat data as 
a function of temperature. Very clear specific heat jump 
is seen at the temperature close to 155 K (see the ex- 
panded plot in the upper inset). A step- like, although 
broad, feature in the specific heat data suggests that 
the anomaly is caused by a second order phase transi- 
tion. A good linear T 2 dependence indicates that the 
specific heat C is mainly contributed by electrons and 
phonons. The fit yields the electronic coefficient 7 =3.7 
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FIG. 1: (a) The electrical resistivity vs temperature for a se- 
ries of LaOi_a;Fa;FeAs. Inset: The phase diagram showing the 
anomaly (red cricle) and superconducting transition (black 
square) temperatures as a function of F content, (b) specific 
heat vs temperature curve. Upper inset: the expanded region 
near the transition temperature; lower inset: a plot of C/T vs 
T 2 . (c) and (d) T-dependent resistivity under magnetic field 
for the pure and 2% F-doping samples, respectively. 
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FIG. 2: (Color online) (a) The reflectance spectra in the 
far-infrared region at different temperatures for the pure 
LaOFeAs sample. The inset shows the reflectance over broad 
frequency range at room temperature, (b) The conductivity 
spectra at different temperatures. 



mJ/mol-K 2 and the Debye temperature #£>=282 K. Note 
that the electronic coefficient 7 is significantly smaller 
than the values obtained from the band structure cal- 
culations being about 5.5-6.5 mJ/rnol-K 2 .[(| 0] This re- 
sult is unconventional, because usually the band struc- 
ture calculation gives smaller value than the experimen- 
tal data, their difference is ascribed to the renormaliza- 
tion effect. We recently measured a Ni-based oxypnictidc 
La(Oo.gFo.i)NiAs, and indeed found that the experimen- 
tal electronic coefficient is larger than the band structure 
calculation [10]. As we shall show below that a partial 
energy gap formation is revealed at the phase transition 
and explained as originated from the SDW instability, 
the smaller experimental value here could be naturally 
accounted for by the gap formation which removes parts 
of the density of states below the phase transition. 

We also measured the resistivity under magnetic field. 
Figure 1 (c) and (d) shows the results for the pure and 
2% F-doping samples, respectively. The anomaly tran- 
sition temperature is rather insensitive to magnetic field 
up to 14 T, however sizable positive magnetoresistance 
are observed at low temperature for pure LaOFeAs. The 
positive magnetoresistance could be understood in terms 



of the suppression of SDW order (and therefore enhanced 
spin scattering) by external magnetic field. After 2% F- 
doping, the low temperature magnetoresistance becomes 
much weaker, which suggests the same tendency that F- 
doping tries to suppress the anomaly around 150K. The 
overall picture is pretty much similar to the elemental Cr, 
a typical SDW system, where transition temperature is 
extremely insensitive to magnetic field and sizable mag- 
netoresistance exists after SDW transition. Q 

Important information about the nature of the phase 
transition could be obtained from the optical measure- 
ment. Figure 2 (a) shows the temperature-dependence 
of the reflectance in the far-infrared region for the pure 
LaOFeAs sample. The reflectance at 300 K over broad 
frequency range is shown in the upper inset. Even though 
the measurement was performed on polycrystalline sam- 
ples, two important experimental findings regarding to 
the phase transition could be unambiguously drawn from 
the measurement. First, below the phase transition tem- 
perature, the reflectance is strongly suppressed below 
the frequency of 400 cm -1 . This is dramatically differ- 
ent from the F-doped LaOo.gFo.i-^FeAs superconducting 
sample where a monotonous increase of the reflectance 
was seenfllj]. The suppression is a strong indication 
for the formation of an energy gap in the density of 
state. Note that the reflectance at very low frequency 
increases fast and exceeds the values at high tempera- 
ture, e.g. above the phase transition. This indicates that 
the compound is still metallic even below the phase tran- 
sition, being consistent with the dc resistivity measure- 
ment which reveals an enhanced conductivity. The data 
indicate clearly that only partial or some of the Fermi 
surfaces are gapped. We note that the gap is very small, 
roughly in the range of 150^350 cm" 1 . Below 50 cm" 1 , 
we find a sharp upturn for conductivity spectrum at 8 K 
(see fig. 2 (b)), suggesting development of very narrow 
Drude compound, which should be linked to the survived 
FS in the SDW state. 

Second, a number of pronounced phonon structures are 
seen in the reflectance and conductivity spectra. Note 
that infrared spectroscopy probes phonon modes only 
near r point, the frequencies of the observed infrared 
active modes are basically in agreement with the cal- 
culations by Singh and Du.[6| Upon cooling the sample 
below the transition temperature, no new phonon mode 
is observed, [jjj Considering the polycrystalline nature of 
the sample which has random orientations of the ab-planc 
and the c-axis, the absence of the any new phonon mode 
strongly suggests that there is no, or at least negligible, 
lattice distortion across the phase transition. Therefore, 
the phase transition could not be of CDW origin. 

In the following we shall illustrate that the transition 
is caused by the SDW transition due to the Fermi surface 
nesting from first-principles calculations. Since only the 
magnetic translational symmetry has been broken, there 
is no modulation for the charge density, which can ex- 
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FIG. 3: (a) and (b) are the Fermi surfaces (FS) cutted at the 
k z —0 and ir planes respectively. The blue lines are original 
calculated FS, while red lines indicate the FS shifted by q=(ir, 
7r, 0). Significant nesting is clear, (c) The calculated Lindhard 
response functionn Xo(o)- It is strongly peaked at M point, 
(d) The Xo(q) along the F-M line. The peak at M point is 
much suppressed by F-doping. 



plain the absence of new phonon modes appearing in the 
optical data after the transition. As we will show below, 
after the SDW transtion, a stripe like spin ordering pat- 
tern appears which breaks the 4-fold rotational symmetry 
and induces a nematic oder in the charge sector. 

Using the experimental high temperature structure 
(P4/nmm symmetry), our calculated electronic struc- 
tures and Fermi surfaces are identical to those done by 
other people @, H 01- B y cutting the BZ into fixed-fc 2 
planes, circle like FS are resolved (as shown in Fig. 3 (a) 
and (b) for k z =0 and k z =ir planes respectively). Surpris- 
ingly, by shifting the circles around the M points to the 
r point, i.e. by a vector q=(ir, ir, 0), the electron-like FS 
will largely overlap with the hole- like FS, suggesting sig- 
nificant nesting effect. For example, for the k z —0 plane, 
the small electron circle almost exactly overlaps with the 
large hole circle after shifting by q. Similar thing happens 
for the k z =it plane. The nesting effect can be quanti- 
tatively estimated by calculating the Lindhard response 
function as shown in Fig 3(c) and (d). The calculated 
Xo(q) is strongly peaked at M point for undoped com- 
pound, and it is much suppress by electron-doping, be- 
cause the up-shift of Fermi level tends to reduce the size 
of hole-like FS and enlarge the electron FS. The exis- 
tence of strong nesting effect would suggest that certain 
kinds of ordering, either CDW or SDW, may develop 
at low temperature. As mentioned above, the absence 
of any new phonon mode after the transition suggests 
that CDW modulation is unlikely, otherwise the struc- 
tural distortion would lead to much more new phonon 
modes. Nevertheless, here we provide further theoretical 
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FIG. 4: (a) The total DOS of (lxl) NM and (a/2 x y/2) AF 
solutions. The reduction of N(Ef) is about 7. (b) The pro- 
jected DOS of Fe site for the AF solution, (c) The stripe-like 
ordering pattern of Fe plane and the charge distribution (oc- 
cupied). The Fe atoms are indicated by red and blue spheres 
with up and down spin respectively. The original (lxl) unit 
cell is indicated by the dashed lines, (d) The calculated FS of 
AF state. The BZ special points of (v2 x v2) cell are indi- 
cated by symbols with prime. The red lines are rotated from 
blue lines, corresponding to the 90 degrees rotation of stripe 
orientation. 



evidences. 

The commensurate nesting vector suggests the dou- 
bling of the unit cell to be (y/2 x y/2) (i.e. the folding of 
the BZ along the dashed lines shown in Fig. 3 (a) and 
(b)), which is not treated by previous calculations 0,0]. 
We have tried to optimized the structure of non-magnetic 
(NM) state without using any symmetry, the solution al- 
ways recover to the original (lxl) structure. However, 
among several possible magnetic solutions, one antifer- 
rromagnetic (AF) ordered state (see the ordering pattern 
shown in Fig. 4) can be stabilized with about 40 meV/Fe 
lowering in energy and about 1.5 ^e/Fe for the spin mo- 
ment compared to (lxl) NM solution. Clearly the nest- 
ing instability likes to couple with spin rather charge. It 
is important to note that only part of the FS are nested, 
as the results, the system remains to be metallic and only 
partial gap opens at Ep. The density of state at Ep is 
reduced significantly compared with original (lxl) NM 
solution as shown in Fig. 4. All those factors are exactly 
what were observed in our experiments. Treating the F- 
doping by virtue crystal approximation, we found that 
both nesting effect and the stabilization energy of AF 
solution are suppressed, again explain the experimental 
observation. 

The obtained AF ordered state is different with the 
(lxl) NM state in the sense that 4-fold rotational sym- 
metry is broken due to the presence of magnetic ordering, 
and the stripe-like ordering pattern is revealed for the Fe 
plane (as shown in Fig. 4 (c)). However, if the stripe 
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orientation is rotated by 90 degrees for the neighboring 
Fe plane stacking along z, the lattice parameters may 
remain to be 4- fold invariant. Furthermore, if we con- 
centrate on the charge (rather than magnetic density) 
distribution, the translational symmetry of total charge 
remains to be the same as original (lxl) structure (as 
shown in Fig. 4(c). As we mentioned, this is the rea- 
son why no new phonon modes are observed. It deserves 
to remark that, in the ordered state, the charge at Fe 
site has a preferentially aligned distribution along stripe 
direction (see Fig. 4 (c)). The fact that charge den- 
sity picks up of preferred direction after a spontaneous 
symmetry-breaking in spin sector could be considered as 
a formation of nematic order. [l3} Here, the order is char- 
acterized by a 4-fold rotational symmetry broken. 
Experiment: 

The samples were prepared by the solid state reaction 
using Fe2C>3, Fe, La, As, and LaF 3 as starting materi- 
als. LaAs was obtained by reacting La chips and As 
pieces at 500 °C for 15 hours and then 850 °C for 2 
hours. The raw materials were thoroughly mixed and 
pressed into pellets. The pellets were wrapped with Ta 
foil and sealed in an evacuate quartz tube under argon 
atmosphere. They were then annealed at 1150 °C for 
50 hours. The phase purity was checked by a powder 
X-ray diffraction method using Cu Ka radiation at room 
temperature. The XRD patterns are well indexed on 
the basis of tetragonal ZrCuSiAs-type structure with the 
space grounp P4/nmm. The electrical resistivity was 
measured by a standard 4-probe method. The specific 
heat measurement was carried out using a thermal re- 
laxation calorimeter. The heat capacity of addenda were 
carefully calibrated before measurement. All these mea- 
surements were preformed down to 1.8K in a Physical 
Property Measurement System(PPMS) of Quantum De- 
sign company. Optical reflectance measurement were 
performed on Bruker 113v and 66v/s spectrometers in 
the frequency range from 20 cmT x to 15,000 cmT x at 
different temperatures. The samples were polished and 
shinny and metallic bright surface was obtained. An in- 
situ gold overcoating technique was used for the exper- 
iment. Optical conductivity was derived from Kramers- 
Kronig transformation of reflectance. The first-principles 
calculations were done using plane- wave pesudopotential 
method and generalized gradient approximation (GGA) 
for the exchange-correlation potential. For the lxl unit 
cell, the calculated results are identical to those per- 
formed by other people @, 0] • 
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